Peak-agnostic high-resolution cis-regulatory circuitry mapping using single cell multiome data.

Single same cell RNAseq/ATACseq multiome data provide unparalleled potential to develop high resolution maps of the cell-type specific transcriptional regulatory circuitry underlying gene expression. We present CREMA, a framework that recovers the full cis-regulatory circuitry by modeling gene expression and chromatin activity in individual cells without peak-calling or cell type labeling constraints. We demonstrate that CREMA overcomes the limitations of existing methods that fail to identify about half of functional regulatory elements which are outside the called chromatin 'peaks'. These circuit sites outside called peaks are shown to be important cell type specific functional regulatory loci, sufficient to distinguish individual cell types. Analysis of mouse pituitary data identifies a Gata2-circuit for the gonadotrope-enriched disease-associated Pcsk1 gene, which is experimentally validated by reduced gonadotrope expression in a gonadotrope conditional Gata2-knockout model. We present a web accessible human immune cell regulatory circuit resource, and provide CREMA as an R package.

Integrated single-cell multiome analysis reveals muscle fiber-type gene regulatory circuitry modulated by endurance exercise.

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.

Steroidogenic Factor 1 Regulates Transcription of the Inhibin B Coreceptor in Pituitary Gonadotrope Cells.

The inhibins control reproduction by suppressing follicle-stimulating hormone synthesis in pituitary gonadotrope cells. The newly discovered inhibin B coreceptor, TGFBR3L, is selectively and highly expressed in gonadotropes in both mice and humans. Here, we describe our initial characterization of mechanisms controlling cell-specific Tgfbr3l/TGFBR3L transcription. We identified two steroidogenic factor 1 (SF-1 or NR5A1) cis-elements in the proximal Tgfbr3l promoter in mice. SF-1 induction of murine Tgfbr3l promoter-reporter activity was inhibited by mutations in one or both sites in heterologous cells. In homologous cells, mutation of these cis-elements or depletion of endogenous SF-1 similarly decreased reporter activity. We observed nearly identical results when using a human TGFBR3L promoter-reporter. The Tgfbr3l gene was tightly compacted and Tgfbr3l mRNA expression was essentially absent in gonadotropes of SF-1 (Nr5a1) conditional knockout mice. During murine embryonic development, Tgfbr3l precedes Nr5a1 expression, though the two transcripts are fully colocalized by embryonic day 18.5 and thereafter. Collectively, these data indicate that SF-1 directly regulates Tgfbr3l/TGFBR3L transcription and is required for postnatal expression of the gene in gonadotropes.

Multi-omics profiling of single nuclei from frozen archived postmortem human pituitary tissue.

Concomitant profiling of transcriptome and chromatin accessibility in isolated nuclei can reveal gene regulatory control mechanisms in health and disease. We report a single nucleus multi-omics analysis protocol optimized for frozen archived postmortem human pituitaries that is also effective for frozen ovine and murine pituitaries and human skeletal muscle biopsies. Its main advantages are that (1) it is not limited to fresh tissue, (2) it avoids tissue dissociation-induced transcriptional changes, and (3) it includes a novel, automated quality control pipeline. For complete details on the use and execution of this protocol, please refer to Ruf-Zamojski et al. (2021) and Zhang et al. (2022).

Transcription factor GATA2 may potentiate follicle-stimulating hormone production in mice via induction of the BMP antagonist gremlin in gonadotrope cells.

Mammalian reproduction depends on the gonadotropins, follicle-stimulating hormone (FSH), and luteinizing hormone, which are secreted by pituitary gonadotrope cells. The zinc-finger transcription factor GATA2 was previously implicated in FSH production in male mice; however, its mechanisms of action and role in females were not determined. To directly address GATA2 function in gonadotropes, we generated and analyzed gonadotrope-specific Gata2 KO mice using the Cre-lox system. We found that while conditional KO (cKO) males exhibited ∼50% reductions in serum FSH levels and pituitary FSHß subunit (Fshb) expression relative to controls, FSH production was apparently normal in cKO females. In addition, RNA-seq analysis of purified gonadotropes from control and cKO males revealed a profound decrease in expression of gremlin (Grem1), a bone morphogenetic protein (BMP) antagonist. We show Grem1 was expressed in gonadotropes, but not other cell lineages, in the adult male mouse pituitary. Furthermore, Gata2, Grem1, and Fshb mRNA levels were significantly higher in the pituitaries of WT males relative to females but decreased in males treated with estradiol and increased following ovariectomy in control but not cKO females. Finally, we found that recombinant gremlin stimulated Fshb expression in pituitary cultures from WT mice. Collectively, the data suggest that GATA2 promotes Grem1 expression in gonadotropes and that the gremlin protein potentiates FSH production. The mechanisms of gremlin action have not yet been established but may involve attenuation of BMP binding to activin type II receptors in gonadotropes, facilitating induction of Fshb transcription by activins or related ligands.

Single nucleus transcriptome and chromatin accessibility of postmortem human pituitaries reveal diverse stem cell regulatory mechanisms.

Despite their importance in tissue homeostasis and renewal, human pituitary stem cells (PSCs) are incompletely characterized. We describe a human single nucleus RNA-seq and ATAC-seq resource from pediatric, adult, and aged postmortem pituitaries (snpituitaryatlas.princeton.edu) and characterize cell-type-specific gene expression and chromatin accessibility programs for all major pituitary cell lineages. We identify uncommitted PSCs, committing progenitor cells, and sex differences. Pseudotime trajectory analysis indicates that early-life PSCs are distinct from the other age groups. Linear modeling of same-cell multiome data identifies regulatory domain accessibility sites and transcription factors that are significantly associated with gene expression in PSCs compared with other cell types and within PSCs. We identify distinct deterministic mechanisms that contribute to heterogeneous marker expression within PSCs. These findings characterize human stem cell lineages and reveal diverse mechanisms regulating key PSC genes and cell type identity.

TGFBR3L is an inhibin B co-receptor that regulates female fertility.

Follicle-stimulating hormone (FSH), a key regulator of ovarian function, is often used in infertility treatment. Gonadal inhibins suppress FSH synthesis by pituitary gonadotrope cells. The TGFß type III receptor, betaglycan, is required for inhibin A suppression of FSH. The inhibin B co-receptor was previously unknown. Here, we report that the gonadotrope-restricted transmembrane protein, TGFBR3L, is the elusive inhibin B co-receptor. TGFBR3L binds inhibin B but not other TGFß family ligands. TGFBR3L knockdown or overexpression abrogates or confers inhibin B activity in cells. Female Tgfbr3l knockout mice exhibit increased FSH levels, ovarian follicle development, and litter sizes. In contrast, female mice lacking both TGFBR3L and betaglycan are infertile. TGFBR3L's function and cell-specific expression make it an attractive new target for the regulation of FSH and fertility.

Single nucleus multi-omics regulatory landscape of the murine pituitary.

To provide a multi-omics resource and investigate transcriptional regulatory mechanisms, we profile the transcriptome, chromatin accessibility, and methylation status of over 70,000 single nuclei (sn) from adult mouse pituitaries. Paired snRNAseq and snATACseq datasets from individual animals highlight a continuum between developmental epigenetically-encoded cell types and transcriptionally-determined transient cell states. Co-accessibility analysis-based identification of a putative Fshb cis-regulatory domain that overlaps the fertility-linked rs11031006 human polymorphism, followed by experimental validation illustrate the use of this resource for hypothesis generation. We also identify transcriptional and chromatin accessibility programs distinguishing each major cell type. Regulons, which are co-regulated gene sets sharing binding sites for a common transcription factor driver, recapitulate cell type clustering. We identify both cell type-specific and sex-specific regulons that are highly correlated with promoter accessibility, but not with methylation state, supporting the centrality of chromatin accessibility in shaping cell-defining transcriptional programs. The sn multi-omics atlas is accessible at snpituitaryatlas.princeton.edu.

Innate Immune Response to Influenza Virus at Single-Cell Resolution in Human Epithelial Cells Revealed Paracrine Induction of Interferon Lambda 1.

Early interactions of influenza A virus (IAV) with respiratory epithelium might determine the outcome of infection. The study of global cellular innate immune responses often masks multiple aspects of the mechanisms by which populations of cells work as organized and heterogeneous systems to defeat virus infection, and how the virus counteracts these systems. In this study, we experimentally dissected the dynamics of IAV and human epithelial respiratory cell interaction during early infection at the single-cell level. We found that the number of viruses infecting a cell (multiplicity of infection [MOI]) influences the magnitude of virus antagonism of the host innate antiviral response. Infections performed at high MOIs resulted in increased viral gene expression per cell and stronger antagonist effect than infections at low MOIs. In addition, single-cell patterns of expression of interferons (IFN) and IFN-stimulated genes (ISGs) provided important insights into the contributions of the infected and bystander cells to the innate immune responses during infection. Specifically, the expression of multiple ISGs was lower in infected than in bystander cells. In contrast with other IFNs, IFN lambda 1 (IFNL1) showed a widespread pattern of expression, suggesting a different cell-to-cell propagation mechanism more reliant on paracrine signaling. Finally, we measured the dynamics of the antiviral response in primary human epithelial cells, which highlighted the importance of early innate immune responses at inhibiting virus spread.IMPORTANCE Influenza A virus (IAV) is a respiratory pathogen of high importance to public health. Annual epidemics of seasonal IAV infections in humans are a significant public health and economic burden. IAV also causes sporadic pandemics, which can have devastating effects. The main target cells for IAV replication are epithelial cells in the respiratory epithelium. The cellular innate immune responses induced in these cells upon infection are critical for defense against the virus, and therefore, it is important to understand the complex interactions between the virus and the host cells. In this study, we investigated the innate immune response to IAV in the respiratory epithelium at the single-cell level, providing a better understanding on how a population of epithelial cells functions as a complex system to orchestrate the response to virus infection and how the virus counteracts this system.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity.

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.